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We investigate neutrino emission from gamma-ray bursts (GRBs) under alternative scenarios 
for prompt emission (the photospheric and synchrotron self-Compton scenarios) rather than the 
classical optically thin synchrotron scenario. In the former scenario, we find that neutrinos from 
the pp reaction can be very important at energies < (10 — 100) TeV. They may be detected by 
IceCube/KM3Net and useful as a probe of baryon acceleration around/below the photosphere. In 
the latter scenario, we may expect ~ EeV p7 neutrinos produced by soft photons. Predicted spectra 
are different from that in the classical scenario, and neutrinos would be useful as one of the clues to 
the nature of GRBs (the jet composition, emission radius, magnetic field and so on) . 

PACS numbers: 98.70.Rz, 95.85.Ry 



I. INTRODUCTION 



Prompt high-energy neutrino emission from gamma- 
ray bursts (GRBs) was predicted in the internal shock 
(IS) model ^ Q and has been studied by many au- 
thors [3, 3, Q (for afterglows, see Ref. [g] and refer- 
ences therein). However, despite recent progresses in the 
Swift era, the mechanism of prompt emission has not 
been well understood Q- One of the most frequently 
discussed scenarios as a standard one is the classical op- 
tically thin synchrotron (including diffusive synchrotron) 
scenario, where prompt photons around the peak energy 
in the hard x-ray or gamma-ray band come from elec- 
trons accelerated at internal shocks and/or by magnetic 
reconnection in the optically thin region. However, this 
scenario cannot satisfactorily explain some observational 
features such as the lower energy spectral index and ob- 
served spectral correlations Q , which may also be related 
to the cooling and efficiency problems 0. Motivated by 
these problems, photospheric emission models have been 
developed [13, fll|l. w hich have an advantage to stabilize 
the peak energy [10, [l^l . 

Another alternative scenario is the synchrotron self- 
Compton (SSC) scenario, which may provide more vi- 
able parameter sets for prompt emission, compared to 
the classical scenario [ij| . The synchrotron peak is in the 
optical/ultraviolet range and gamma-ray photons arise 
from inverse Compton scatterings. This scenario may 
explain some bursts such as GRB 080319B [14], despite 
the lack of bright optical emission in many bursts [l5| . 

In this work, we investigate neutrino emission under 
photospheric and SSC scenarios. Our method of calcu- 
lation based on Geant4 is basically the same as in Refs. 
[3|, |6| , but improved qualitatively and quantitatively by 
including the pp reaction and all the relevant processes 
of protons, mesons and muons [1^1 [l3- ^^ for GRB neu- 
trinos, cooling of mesons and muons is remarkably im- 
portant, which makes neutrino spectra complicated and 
affects estimate of event rates 0, 0, [l3 ■ Especially, we 
demonstrate that pp neutrinos can be important below 



100 TeV. 



II. THE PHOTOSPHERIC SCENARIO 

GRB prompt emission is considered to be radiation 
from a highly relativistic jet toward us, and the typical 
Lorentz factor of the collimated outflow is F ^ 10^^''. 
In the standard picture, a significant fraction of the out- 
flow energy is converted to the radiation energy via some 
dissipation mechanism within the outflow (e.g., inter- 
nal shocks or magnetic reconnection), leading to the ob- 
served isotropic photon luminosity L^ ^ 10^^ ergs s~^ 
[7|. In photospheric emission models, the prompt emis- 
sion comes from around the photospheric radius r ~ Tph, 
at which the Thomson optical depth T-ye — ?^eCT(?'/r) 
is unity. The photospheric radius is given by Tph ~ 



(LA,/crT/47rr^TOpC^) 



1.2 X 10^^cmLM,52.5r2.5, when 
the outflow is baryonic (rig ~ rip). Here Lm is the 
isotropic outflow luminosity carried by baryons, in the 
observer frame 0, [l^]. Possibly, the outflow may be 
dominated by pairs. If rie/np '^ mp/rrie, we obtain 
reh « {LM(TT/4.TrT^meC^) ~ 2.2 x 10^^ cmLAf,52.5r^5 
[ll|. In this scenario, dissipation and thermalization 
should occur around/below the photosphere. For in- 
stance, in the dissipative photospheric model [ll|, we can 
assume that the dissipation is maintained after the coast- 
ing radius, leading to the temperature /cTob ~ 100 keV in 
the observer frame. The observed peak energy e^,^ ^ 500 
keV can be achieved by the Comptonized thermal pho- 
tons [l3|. 

As an example, let us consider an internal collision of 
two subshells in the relativistic outflow, with Tg ^ 10^ 
and Tf ~ 10'^. The Lorentz factor of the merged subshell 
is estimated as F « ^y^f^s — 10^'^, and the Lorentz fac- 
tor of the internal shock in their center of mass frame 
is - V^ « V(r//r. + F,/F/)/2 ~ a few. Hence, 
internal shocks would be mildly relativistic shocks, at 
which electrons can be accelerated, and prompt emission 
mainly comes from the shocked subshells which could 



be magnetized via plasma instabilities at collisions Q- 
Not only electrons but also baryons can be accelerated 
even around/below the photosphere [3, [l3]- The ac- 
celeration time scale of protons in the comoving frame 
of the merged subshcUs is written as iacc — rjSp/ eBc. 
Although realistic values of 77 are not well known, we 
could expect ry ^ (1 — 10) in the most efficient cases, 
when we consider the Fermi acceleration mechanism in 
the Bohm limit 0, [iMl- Here we assume 77 ~ 10 as 
in Refs. [3, Q. The magnetic field is written as B « 
8.4 X IQ^Geii^Lll^^^T-lr-l^^,, where ^b ^ Ub/U^- A1- 
though ^B is uncertain, we expect that the equiparti- 
tion between the magnetic energy density Ub and pho- 
ton energy density U-y can be achieved in the photo- 
spheric scenario, as often expected in the classical one. 
The maximum energy of protons is set by the condition 
ta.cc < niin[ip, idyn]- Here tdyn ~ r/Tc is the dynamical 



hjl^{TeB/m,c) ~ 20 eYel^tJUL. 



1/2 ^1/2 



time scale and t„ ^ = iRH+O+O+^svn+^iC +^=,i+i, 



1 

esc 



is the proton loss time scale (in the comoving frame), 
where tsH, ip7, tpp, tsyn, he, iad, and icsc are coofing 
times of the Bethe-Heitler process, photomeson produc- 
tion, pp reaction, synchrotron emission, inverse Compton 
emission, adiabatic expansion, and the escape time in 
the Bohm diffusion approximation [3, i9|. Following Refs. 
[Sj, 6], we evaluate the maximum proton energy for each 
parameter set, whose typical values are E^'^^ ~ 10^"^*^ 
GeV, in the observer frame. 

High-energy protons can interact with photons and 
protons, producing mesons and muons that decay to neu- 
trinos. As in Refs. @, S 0, 0i i* is often convenient to 
use the effective photomeson production optical depth, 
fp-y ~ tdyn/tp^y. By using the A-resonance approxima- 
50PeVri5/e^,3i6k,v Siii, we have 

(E'p 



tion at Ep 



fi 



pi\ 



23 



7 '^''"'^u — T^e, where L^(< L^) is the 



observed photon luminosity at e^,^. It implies that almost 
all the energy of high-energy protons can be used for pho- 
tomeson production when tp^ is the most important, and 
accelerated protons will be depleted. Similarly, the effec- 
tive optical depth for pp reaction is fpp « KppUpappl ~ 
0.05r-ye, where app ~ 5 x 10~^^ cm~^ and Kpp ~ 0.5 — 0.6. 
The meson production efficiency /meson can be approx- 
imated by ~ min[l,max(/p-y, /pp)], as long as tp^ or tpp 
is the most relevant loss time scale. We can expect that 
a significant fraction of proton energy is used for me- 
son production, which is demonstrated through numeri- 
cal calculations. 



III. THE SSC SCENARIO 

SSC emission models have often been discussed |13| . 
where the observed peak energy e^j^ ~ 500 keV is identi- 
fied with the second peak formed by up-scattered syn- 
chrotron photons. Rather large emission radii close 
to the deceleration radius, r^ec "^ 10^^ cm (^ rph ^ 
10^^ cm) are expected [Ij, |lj|. The first synchrotron 
peak in the observer frame is estimated as e^jj^ « 



B, -2-^7, 52' 15.5' 



and then 



the second peak coming from inverse Compton scatter- 



ings is e^l 



7^ s''i 



700 keYei_,eB%L\%rY,\. In 



addition, the third peak of e^J^ ~ 20 GeV can also be 
expected, which is one of the predictions in this scenario. 
Most of the radiation energy can be released as high- 
energy gamma rays, so that we expect £l^° ~ Y£^2 with 

y ~ ^g , where Y is the Compton Y parameter [?[ and 
£1^2 is the isotropic energy around the second peak. 

Note that, in this scenario, the magnetic field should 
be smaller than in the equipartition case. Correspond- 

10-2 leads to B - 



g _2-^j 52^ 2 5''i5"5- Then, from the condition tacc < 



ing to r ~ 10, ^B = Ub/U^ 

min[ip, idyn], the maximum proton energy is typically es- 
timated as S™^'' ~ 10^"-" GeV for rj ~ 10. Production 
of ultra-high-energy cosmic rays (UHECRs) may be pos- 
sible. 

We can also evaluate the effective photomeson optical 

depth as fpjiEp) ~ ^7-77^ 

'15-5^ 2.a-oD,ai!cv 

below £;" 



at sufficiently high energies ^v.^^^.. ^p 

500EeVr2.5/£ob32cv When tpj is the most rele- 
vant, we expect a significant fraction of proton energy is 
used for the photomeson production around the highest 
energies. But, in this scenario, fpj and fpp are rather 
small at lower energies, so that efficient < PeV neutrino 
production is not expected. 



IV. THE NEUTRINO SPECTRUM AND FLUX 

In this work, we numerically calculate neutrino spec- 
tra and evaluate fluxes through the method used in Refs. 
[3, [y]. Calculations are first executed in the comoving 
frame, and then results in the observer frame can be eas- 
ily obtained via the transformation. The photomeson 
production and pp reaction are treated in detail, as well 
as various cooling processes of mesons and muons, i.e., 
synchrotron, inverse Compton, meson-photon, TTp/fj,p 
and adiabatic cooling processes. Neutrino spectra can 
be calculated, giving a target photon spectrum with 
Uj — (Lj/^TTr'^T'^c), target (thermal) proton density 
Hp — {Ut\i/mpC^) and Ub = £,bUj- For the evaluation 
of neutrino fluxes, proton fluxes should also be given and 
we use dncR/dSp oc e" ^ in this work. Although the pro- 
ton spectral index of p ~ 2 is often expected for nonrela- 
tivistic or mildly relativistic shocks with the compression 
ratio of ~ 4, different values are possible, for example, 
due to large angle scatterings rather than small pitch- 
angle scatterings across relativistic shocks [la, [3 • The 
baryon energy input fcR is given by using the nonther- 
mal baryon loading factor ^acc = ^ckI^-j = '^cr/^7°' 
where £^ and £^^° are the geometrically corrected radia- 
tion energy and isotropic radiation energy respectively. 
Here we adopt ^acc ^ 1 — 10, as in Refs. QjHHilal- Differ- 
ences in spectra mainly come from r, S,b and the photon 
spectrum, which depend on prompt emission scenarios. 
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FIG. 1: The (Vfj.+i'fj,) fluence from a GRB event at z — 0.1. pp 
neutrinos (thick line) and py neutrinos (thin line) are shown. 
Baryonic photosphere: r — rph = lO'^^'^ cm {r-ye = 1); Nf^ ~ 
1.7 events. Baryonic subphotosphere: r — 10~'''^rph (tje = 
10); A^^ ~ 1.4 events. Pair photosphere: r = rph ~ 10^*'^ 
cm {r-ys = 1); A'^^^ ~ 0.018 events. Pair subphotosphere: r — 
10"°-^reh (r^e = 10); A^^ ~ 0.16 events. 4'° = lO^^'^ ergs, 



fph 



' 'ye 



r — 10^'^, ^B ~ i, and ^acc ~ 1 are adopted. We assume 
Tie — Up in case of the baryonic photosphere while rie — ^^^Up 
in possible cases of the pair photosphere. 



For instance, in the classical scenario, r ^ 10^'^ ^'"^ cm, 
^s ~ 1 and the broken power-law photon spectrum are 
typically used 0, i, i, [1] . 

In the photospheric scenario, we may expect smaller 
radii of r ^ 10^^^^^ cm and strong magnetic fields with 
^_B ~ 1. Prompt emission comes from the photosphere 
where r-ye = 1, so that we can adopt the broken power- 



law photon spectrum as dn/de oc {s/e^) for e" 



< e < 



/3 



e and dn/de ex (e/e ) for e < e < s"^^'^, expressed 
in the comoving frame. Here we set e™'" = 1 eV and 
ginax _ 2^ MeV, because the synchrotron self-absorption 
and pair-creation absorption will be crucial below and 
above these energies, respectively. But our results are not 
sensitive to them. We take e^ — s'^^^/T = 1 keV, a — 1 
and f3 = 2.2, which are obtained from observations. How- 
ever, since photons would be significantly thermalized 
when Tje > 1, we adopt the black-body photon spectrum 
with the temperature T — [U-^/ a) for subphotospheric 
emission. Up is given assuming C/th = U.y , and we use r^e 
as a parameter instead of L^ . 

In the SSC scenario, we expect larger radii of r ^ 
^gi5-i6 ^j^ ^^^ relatively weak magnetic fields with 
^B ^ 0.01. The photon spectrum is given by the sum 
of multi-broken power-law spectra with break energies, 
e" = (10-2 - 10-1) eV and e^'^ = 1 keV with a = 1 and 
(3 — 2.2. Note that larger emission radii lead to lower syn- 
chrotron self-absorption energies of e™'" = (10-^ — 10- 1) 
eV. For comparison, we also show neutrinos produced via 
the p7 reaction between protons accelerated at the ex- 
ternal reverse-shock (RS) and prompt photons produced 



FIG. 2: As in Fig. 1, but z = 0.9. SSC (IS-IS) (thick 
solid line): r = 10^^ cm, V = 10^ U 



7 



10'^^■^ergss-^ 



2 X 10 

-,1B.4 



■1 events. SSC (RS-IS) (thin 
cm, Px ~ 150, and BC ~ 7.7 G 



and e,B = 0.01; A^, 
solid line): rx — 10 
(coming from E]^° = 4 x lO^** ergs, Aq = 7.5 x lO" cm, 
^b/Ib = 4 X 10-^ and 4. = 0.4, inferred in Refs. [ll|); 
N^, ~ 10-^ events. For SSC, 42° = 10*''"^ ergs, a = 0.86, 
13 = 3.6, and ^acc = 1. Photosphere (baryonic) (thick dashed 
line): r = rph = 10^^'^ cm; A'^ ^^ 0.2 events. Photosphere 
(pair) (thin dashed line): r = rph = 10^*'^ cm; A^^ ~ 2 x 10-^ 
events. Subphotosphere (baryonic) (thick dotted line): r = 
lO-^'^rph; A'^ ~ 0.2 events. Subphotosphere (pair) (thin dot- 
ted Hue): r — lO-^'^rph; A'^ ~ 10-^ events. For photosphere 
and subphotosphere, £^^° = 10^^'^ ergs, ^b ~ 1, and ^acc ~ 1. 



by internal dissipation. The relevant quantities at the 
crossing time, such as the radius rx , Lorenz factor of the 
shocked ejecta Fx and magnetic field i?x, are evaluated 
according to the standard reverse-forward shock theory 
7] . The procedure is described in Ref . [G] . 

The results in the photospheric scenario are shown in 
Fig. 1. We can see a pp neutrino component is dominant 
at En < 10 TeV, and more important for subphotospheric 
emission. The Bethe-Heitler and photomeson production 
processes become more relevant at higher energies. The 
former can lead to a dip between pp and P7 components 
(dotted-dashed lines). The latter makes P7 neutrinos, 
whose fluence is suppressed at sufficiently high energies 
since mesons and muons cool before they decay. For ex- 
ample, we can find break energies around PeV, and p^ 
neutrinos from kaons at higher energies E^ >^ 10 PeV, 
for baryonic photosphere and baryonic subphotosphere. 
When the subphotospheric emission is expected, dissipa- 
tion would continue from the subphotospheres all the way 
to the photosphere, enabling us to expect that pp neutri- 
nos from subphotospheres are more important at lower 
energies, while P7 neutrinos from around/above the pho- 
tosphere at higher energies. Muon event rates expected 
by IceCube, A^^(> lO^'^ GeV), are also shown in the fig- 
ure caption. If a burst occurs at 2; < 0.1, we may detect 
a few events. Event rates of pp neutrinos dominate over 
those of p7 neutrinos for baryonic subphotosphere. 



In Fig. 2, we show neutrino spectra in the SSC and 
photospheric scenarios for an energetic burst such as 
GRB 080319B. In the SSC scenario, expected event rates 
are few due to the small photon density at large radii. 
Note that neutrinos coming from protons accelerated at 
the reverse shock are dominant at lower energies [see SSC 
(RS-IS)] . It is because reverse-shock protons interact with 
blue-shifted prompt photons, which are assumed to have 
the steep photon spectrum with /3 = 3.6 (solid lines) p^ . 
In the photospheric scenario, we can expect higher neu- 
trino fluences than in the other scenarios. 

It is important to consider the cumulative neutrino 
background, since the time- and space-coincidence is ex- 
pected for GRB neutrino emission [41. The background 
flux can be roughly estimated as [3, S 01 
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where /sup is the suppression factor due to cooling of 
mesons and muons 0, M, H [13 and fz expresses the 
contribution from the high redshift GRBs [2, Q • Here 
-Rgrb(O) is the overall (not apparent) local rate, where 



m. 



GRBs are regarded as jets with S^ ^ 10 ergs 

The results are shown in Fig. 3, where the GRB 
rate evolution is properly considered [6|. Predictions 
in the photospheric and SSC scenarios are rather differ- 
ent from the classical one. In the photospheric scenario, 
we expect larger event rates than others with the same 
baryon loading, and pp neutrinos may become important 
at < (10 - 100) TeV. In the SSC scenario, we could ex- 
pect '^ (1 — 10) EeV neutrinos from interactions between 
UHECRs and soft photons, but note that our evaluation 
would be applied to only a fraction of GRBs [la] • 



FIG. 3: The cumulative neutrino backgrounds from GRBs. 
Classical: originally predicted by Waxman and Bahcall [g] 
and the numerically calculated spectrum is taken from Refs. 
[3,[g|, but averaged over sets A and B, F = 10^'^ and ^s ~ 1; 
A^^ ~ 9.4 events/yr. Photospheric (baryonic) (thick dashed 
line): r = rph = 10^^'^ cm, F = 10^■^ and ^b = 1; A"^ ~ 



Photospheric (baryonic) (thin dashed line): 



Nu 



16 events/yr. 

-|2.5 



24 events/yr. 

r = 10-°Vph, F = 10^", and ^s = F ^y^ 
SSC (IS-IS) (thick dotted hne): r = lO^^'^ cm, F = 10^ 
L!; = 10'^^-^ergss-\ and ^s = 0.01; A,, ~ 0.14 events/yr. 
SSC (RS-IS) (thin dotted hne): rx ^ lO'^'^ cm, Fx ~ 160, 
and Bx ^ 5.9 G; A^ ~ 0.014 events/yr. Flare/late prompt: 
proposed by Murase and Nagataki [3] and the numerically 
calculated spectrum is taken from Refs. [3, [g| (the model 
LPO); Ap ~ 1.2 events/yr. WB: Waxman-Bahcall bounds 
shown as benchmarks 2]. The cosmic- ray proton energy in- 
put per logarithmic interval iShecr = £cr/R is normalized 
to 10^^ ergs for the prompt emission scenarios while 10^° ergs 
for flare/late prompt. As a GRB rate evolution model, the 
GRBS model in Ref. [g is adopted with the cosmological pa- 
rameters (fim = 0.3, f^A = 0.7; Ho = 71kms"^Mpc"^) and 



V. IMPLICATIONS AND DISCUSSIONS 

In the near future, high-energy neutrino signals may 
be detected by IccCube and KM3Net E^. IceCube Deep 
Core for ~ TeV neutrinos and acoustic/radio and shower 
detectors for ~ EeV neutrinos will also be useful. In 
this work, we have first demonstrated that pp neutri- 
nos can become more important than p"f neutrinos in 
prompt neutrino emission, if the photospheric scenario is 
realized. Detection of < (10 — 100) TeV neutrinos is im- 
portant as a probe of dissipation and baryon acceleration 
around/below the photosphere as well as a diagnosis of 
the jet composition, although the detectability depends 
on the pair loading. Also, these photospheric neutrinos 
can be distinguished from precursor ones [17] through 
correlations of neutrinos with prompt photons. At higher 
energies, pj neutrinos become more important, and there 
may be contributions from optically thin internal shocks 
occurring above the photosphere. However, pj neutrinos 
cannot be expected for too large values of rj, although we 
have assumed 77 ~ 10 so far. When 77 > 10'^"^, only pp 



neutrinos may be relevant. Note that our predicted fluxes 
are below the current AMANDA limit and Waxman- 
Bahcall bounds [2l'| . If neutrinos are not detected in the 
future, it implies that baryon acceleration is insufficient 
(small ^acc and/or very large 77) or that prompt emission 
occurs sufficiently above the photosphere. In addition, 
we have demonstrated that synchrotron photons in the 
optical/ultraviolet range may enhance ~ EeV neutrinos 
in the SSC scenario, compared to the other scenarios. 
Cooling of meson and muons is not so important, and ^ 
EeV neutrinos will be useful as a probe of UHECR ac- 
celeration. Detection of them will also imply the nature 
of GRBs, e.g., r and ^5. 

One may expect hadronic high-energy gamma rays can 
be detected by the recently launched Fermi satellite. For 
photospheric emission, we do not expect that high-energy 
gamma rays with £ob ^ a few xFrngC^ escape from the 
source due to the large optical depth for pair creation 
[22| |. Primary hadronic gamma rays induce the cascade 
in the source, and resulting spectra would be similar to 



those in the classical scenario. However, as indicated 
in Asano and Inoue [23], it is not easy to find observa- 
tional signatures of hadronic gamma rays with ^s ~ 1, 
^acc ~ 1, a = 1 and /3 = 2.2, considered by us. (Note 
that the steeper electron spectral index of pe ^ 3.0 and 
softer low-energy photon index of a ^ 1.5 in their calcu- 
lations could lead to overestimating hadronic signatures 
in spectra [2^.) Furthermore, we would not expect to 
distinguish between the classical and photospheric sce- 
narios from signatures of cascaded hadronic gamma rays 
themselves. Hence, detection of < (10 — 100) TeV pp 
neutrinos is more important as a unique probe. On the 
other hand, in the SSC scenario, primary high-energy 
gamma rays produced via the photomeson production 
could escape from the source. Their typical energy is 
> EeV, and they cannot avoid attenuation by the cos- 
mic background photons. The detectability of cascaded 
gamma rays in the GeV-TeV range strongly depends on 
the intergalactic magnetic field strength. Only when the 
intergalactic magnetic field is very weak. Big ^ 10^^^ G, 
we may detect them as a pair echo 25[. However, even in 
such cases, we expect that the pair echo emission would 
be subdominant compared to the primary leptonic emis- 
sion. In fact, in the SSC scenario, strong high-energy 



gamma rays should be generated through the SSC emis- 
sion, and the isotropic radiation energy around the third 



peak, £-° 



10^ 



54.5 



ergs, is larger than the isotropic en- 



ergy of primary hadronic very high-energy gamma rays, 
'^7 VHE ~ 10^^-^ ergs. Hence, it would be difhcult to use 
hadronic gamma rays as a probe of baryon acceleration 
in the SSC scenario. 

GRBs may be the origin of observed UHECRs, as is 
discussed in many papers. Note that, if the photospheric 
scenario is real, UHECR explanation is impossible due to 
strong proton losses at the highest energies. In the SSC 
scenario, UHECRs may be explained. 
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